The Pseudogap Challenge: Understanding the ab Plane AC Conductivity below T c 
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We establish that the cuprate pseudogap plays a crucial role in the afc-plane optical conductivity 
a(uj,T) for temperatures T < T c . The pseudogap signatures in a(ui,T) associated with competing 
proposals of a "hidden order" and "superconducting" origin for the d x 2_ y 2 pseudogap are found to 
differ in their qualitative U),T dependences. For the latter case, as T increases from 0, excess low ui 
weight appears; moreover, a much wider range of U) contributes to form the condensate. We discuss 
these theories in light of current experiments. 



While a major fraction of the cuprate phase diagram 
has been found to be associated with a pseudogap there 
is still no definitive explanation for the origin of this 
anomalous phase. Here normal state thermodynamic and 
transport properties exhibit an excitation gap of rf-wave 
symmetry (with onset at temperature T*) which appears 
to be continuous with the excitation gap of the super- 
conducting phase (with onset at T c < T*). Theoreti- 
cal schools for addressing this phenomenon fall roughly 
into two camps in which this pseudogap derives from 
(i.e., is "intrinsic") or is independent of (i.e., "extrinsic" 
to) the superconductivity. Despite its clear presence 
in thermodyjiamical and transport data, mysteriously, 
experiments^ on the in-plane electrodynamics show few 
conspicuous signatures of the pseudogap. In this pa- 
per we demonstrate how, in both theoretical schools, the 
pseudogap enters as a rearrangement, relative to BCS 
predictions, of spectral weight in Reer(cj) = <Ti(u>). More- 
over, we suggest that it is the width of the range of fre- 
quencies contributing to form the condensate which pro- 
vides a most effective criterion, both for descriminating 
between theoretical approaches, and for direct compari- 
son with the data. The experimental observationHa that 
(in the underdoped regime), relatively high lu participate 
in the condensate poses a serious challenge for both strict 
BCS theory and the extrinsic school. 

Strict BCS theory is overly constrained because it con- 
tains a single energy scale. Nevertheless, we argue here 
that a mean field theoretic approach is appropriate be- 
cause of (1) the notable similarity of many (but not 
all) measured properties to their BCS predictions. (2) 
the very narrow critical regime, and (3) the large sep- 
aration between T* and T c (which appears to be too 
big to derive from fluctuations around strict BCS the- 
ory). Moreover, the electromagnetic response of the su- 
perconducting state can be readily formulated at this 
generalized mean field level. In the standard formu- 
lation, the imaginary axis conductivity is written as 

<Jab{Q) = [Pab{Q) - n a b] where Q = (q = Q,iu) n ). 

We define Green's functions 



K = (e k - li) 2 + (A k ) 2 



(2) 



Here ek refers to a generalized "band structure" , which 
may differ from the bare hopping dispersion £k = 
— 2t(cosk x + cosky), whereas Ak refers to a general- 
ized "gap" parameter, whose k-dependence we separate 
as Ak = A(^k with <^k = cos(fca;) — cos(fc y ) for d-wave 
pairing. 

The quasiparticle velocity v a — d a e k leads to 



i ab = 2Y,[d a e k d b e k (FF - GG) 
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BCS theory predicts the response Pab CS (Q) — 
-2 Y.k d a e lt d b e k (GG++FF+) (where G = G(K), G+ = 
G(K + Q), and similarly for F,) so that the general ex- 
pression for the superfluid density (n s ) ab = n a b — Pab(0) 
reduces in the BCS case (omitting tensor indices) to 
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Since in BCS theory there is only one energy scale, one 
identifies A as A sc = A, and n s vanishes only when 
the excitation gap vanishes. More correctly, it should be 
clear that the superfluid density must necessarily reflect 
the superconducting order parameter A sc which may be 
distinguished from the total excitation gap A, so that 
n s (T) cx (A* C (T)) 2 . This serves to underline a major 
shortcoming).! of strict BCS theory. 

The simplest mean field pseudogap approach (case (i)) 
is referred to as "rescaled BCS theory" u. Here, the sin- 
gle energy scale framework of BCS theory A sc = A is 
preserved while modifying the response by a frequency- 
independent factor as P a b — l(T)Pab CS - This enhance- 
ment causes n s to vanish at T c lower than T*. The 
factor 7(T) is related to the suppression of the super- 
fluid n s = f(T)n® cs by the sum rule constraint j(T) — 

n-nft-W 1 " /(T)) + /(T) ' With /(T) ChOS6n t0 Van " 

ish at T c and fitted to n s {T). The microscopic physics 
of this otherwise phenomenological scheme is presumably 
connected to the approach of Lee and WenH; 7 may also 
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reflect the renormalization of electron charge as can occur 
in Fegmi-liquid based approaches to the superconducting 
stated. In overdoped samples with T c — T* , BCS theory 
works well, i.e., f(T) = 1. 

Alternatively, one may argue that A sc 7^ A, so that 
at the mean field level, the gap and the order parameter 
arc distinct. The extrinsic school introduces a pseudo- 
gap A pg directly into the quasi-particle bandstructure ek 
and defines A = A sc in the above mean field formal- 
ism, yielding n s oc A 2 C via Eq. (Q). Below T c , Eq. (||) 
shows that in this school the dispersion of the fermionic 
quasi-particles has no simply defined BCS-like excitation 
gap. However, in this paper (as elsewherea) we will as- 
sume that one can define a meaningful excitation gap for 
the extrinsic case as Ak w ^/(A^ 7 ) 2 + (A|. c ) 2 . We im- 
plement this picture by introducing an additional order 
parameter which breaks translation symmetry, leading to 
two bands — ±\/£k + (^k 9 ) 2 - The precise superfluid 
density expression is slightljj-imore complicated than Eq. 
(||) due to interbank effects.^ While we calculate P a b for 
a d-density wavaljlld (DDW) here, we argue that our 
results are generic, applying to other extrinsic pseudo- 
gap theories BEj The Green's functions F^^G^ for the 
two bands are obtained from Eq. (|l|) , (||) by replacing ek 
with e k . For the ei-density- wave/superconducting coex- 
istent state, we find 

P ab {Q) = -2 £ \vPv£\G»G» + + F°'FX) (5) 

K,v=± 

+ v^vl 2) (G v G+ v + F V F+ V ) 

where momentum is summed over the half Brillouin zone, 
v is the band index, = <9 a ek, and = (^d a A^ — 
A k "7A,i k ) <u- 

Finally, we consider a third mean field alternative (iii) 
in which the pseudogap is intrinsically connected with 
the superconductivitytJ. In this picture A pg is associ- 
ated with a stronger-than-BCS attractive interaction, so 
that (finite momentum) pairs form at a higher temper- 
ature T* than that at which they Bose condense (T c ). 
Here the pseudogap enters on the same footing as the 
superconducting order parameter, and the fermionic dis- 
persion is given by Eq. (||) with the value of Ak assigned 
to be Ak, and ek = £k- In ti4s last approach, which has 
a detailed microscopic basisEj, the superfluid density has 
the form 



n s = [A 2 c /A 2 ]t^ cs ; A = \/ (A») 2 + (A sc ) 2 . (6) 

The corresponding electromagnetic response is the sum 
of fermionic (BCS, with full gap A(T)) and bosonic con- 
tributions 

P ab (Q) = -2Y,da&d b UFF++GG+) 

K 

+ A^2d a n p d b n p tP3(p)tps(p + Q) (7) 
p 



Here we have approximatedE-3 the second bosonic term, 
which originates from an Aslamzov-Larkin diagrams 
It represents the coupling of radiation to uncondensed 
pairs of charge-2e responsible for suppressing n s rela- 
tive to nf cs . The pairs have a propagator t pg (P), 
which at small P can be parameterized as t~ g (p,u>) = 
lu(1 + iv) — tt p (after continuation to real frequency 
lo.) Here il p = p 2 /2M and the pair chemical potesu. 
tial is zero for T < T c . While our microsopic theoryts 
yields the pair mass M and lifetime parameter associated 
with z/, the diagrammatic calculations can bs-avoided by 
borrowing from general TDGL formulationsEZI which de- 
scribe the hydrodynamics of non-condensed bosons. In 
particular, the optical sum rule constrains the conduc- 
tivity spectral weight of this boson contribution to be 
n boson _ [A 2 s /A 2 ]n^ cs , determining the prefactor A. 

In the intrinsic school the condensate (A sc ) grows 
with decreasing T at the expense of A pg as the frac- 
tion of (quasi-ideal Bose gas) finite momentum pair ex- 
citations decreases. Precisely at T = 0, A pg vanishes 
and the ground state is given by the BCS generaliza- 
tion first proposed in Ref. [l8|. By contrast, in the ex- 
trinsic model A pg is relatively T independent below T c 
in the underdoped regime. One can sunuaarize: to a 
good approximation, for the intrinsic modell3 (case (hi)) 
A 2 g (T) « A 2 (T C )(T/T C ) 3 / 2 , whereas for the extrinsic 
model (case (ii)) at low doping A 2 g (T) w A 2 g (T c ). 

The following calculations for ui{oj) address an un- 
derdoped (UD) and a slightly over-doped (OD) cuprate 
with hole concentration x quantified by excitation gap 
ratio A(T c )/A(0) which assumes the values 0.99 and 0.5 
respectively. This corresponds roughly to x ~ 0.06 and 
x « .17. The two quantities n s (0)/A(0) and T c /T* were 
computed within each microscopic theory and roughly 
fit to experiment. We presume that the temperature de- 
pendence of the full excitation gap for T < T c is given 
by the BCS d-wave result, so that A(0)/T* is a univer- 
sal number, as seems to be the case experimentallyEl 
We include scattering from random non-magnetic impit. 
rities by renormalizing F and G in the above formulae,! 2 !!! 
and, for simplicity, compute the bosonic term in Eq(7) 
using standard TDGL resultalZl. Following conventional 
procedures,! 2 -^ we pre-calculate a table of self-energies on 
the real frequency axis. We handle the complicated be- 
havior of energy and k integrands using automatic Id 
or 2d adaptive quadrature, determining o\ to 1% rela- 
tive accuracy. While we have explored both Born and 
unitary limits, the computations are more tractible in 
the unitary limit, which we consider here with scatter- 
ing potential u5(r), u/2t — > 00, and with 0.01 impurities 
per unit cell. Our numerical calculations (1) conserve 
the integrated density of states, conserve conductiv- 
ity spectral weight, and (3) yield£3 ui(T = 0, lj = 0) = 
(2/7r)2f/A. Here and throughout we use conductance 
units e 2 /h = 1. 

In the insets of Figures la and 2a are plotted the nor- 
malized superfluid densities for the intrinsic (iii) and ex- 
trinsic cases (ii) in the underdoped regime. The strict 
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FIG. 1: Real part of ac conductivity for the intrinsic (solid), 
rescaled BCS (dotted) and strict BCS (dashed) models at 
underdoped (UD) and slightly overdoped (OD) x at T — 
(la, lc) and T = T c (lb, Id). Insets plot the corresponding 
n s compared to nj cs (la) and the bosonic contribution to o\ 
(lb). 
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FIG. 2: <Ti(w) for the extrinsic psudogap at underdoped (UD) 
and slightly overdoped (OD) x at T = (2a,2c) and T = T c 
(2b,2d). The arrow indicates A(0) « A P9 (T C ) 



"BCS" curves for each case are chosen to have the same 
excitation gap parameter, so that nf cs vanishes at T*. 
It is clear from the above discussion that the calculated 
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FIG. 3: Partial integral of condensing spectral weight, nor- 
malized by the total, for the three pseudogap models (under- 
doped). 



superfluid density is depressed and vanishes at a lower 
temperature (T c ) than its BCS counterpart. Sum rule 
considerations thus indicate that there must be excess 
spectral weight in (T\(uj, T) at cu ^ associated with pseu- 
dogap effects. 

The frequency distribution of this excess spectral 
weight can be inferred from the plots in Figures 1 and 
2 which show <j\ (us) for the UD and OD hole concentra- 
tions at T = and T = T c in the intrinsic (solid lines) 
and extrinsic cases respectively. In Figures la, and lc 
the rescaled BCS model (dotted lines) coincides with the 
solid curves at T = due to the complete condensation 
of the pseudogap. It is clear that in the intrinsic case 
the added pseudogap weight (jajesent at any T ^ 0) ap- 
pears at very low oj, varying asO 1 /y/tJ, as well as at very 
high uj where it introduces a long tail associated with the 
fermion induced damping of non-condensed pairs. Inter- 
actions with the fermionic two particle continuum (which 
itself is gapped) constrains the bosonic lifetime. By con- 
trast, in the extrinsic case the pseudogap is associated 
with a secondary maximum which appears at the en- 
ergy scale A pg (T c ) « A(0). This arises from interband 
effectslill with threshold at lu = 2/i. 

It should be noted that a direct comparison of o~\ (u) 
between theory and experiment is complicated because 
the origin Qf-,the so-called mid-infrared spectral weight 
is uncertain^, and because the ratio of the T = con- 
densate to the total integrated spectral weight (at T c ) 
is anomalously small (20% or less)E3. Thus, we focus 
on the difference curves, and their associated integrals. 
We introduce the partially integrated spectral weight 
AiV(O) = | / f2 [fJi(w,T c ) - a!(u),0)]du), which is plot- 
ted in Figure 3 for the underdoped regime. This plot 
indicates which frequencies are required to account for 
the superfluid density in both the extrinsic and intrinsic 
cases. 

Strikingly, the intrinsic pseudogap and rescaled BCS 
theories both show substantial contributions to the con- 
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densate up to relatively high frequencies of order 4 
A(0) « 960cm _1 or more. By contrast, in the extrin- 
sic school, the condensate derives from a much narrower 
range of frequencies. Since intcrband conductivity terms 
make a negligible condensate contribution, this observa- 
tion holds evieja for extrinsic theories which lack inter- 
band effects@ll3. An analogous narrowing of condensate 
frequencies occurs in strict BCS theory at low x when 
A(T) is properly chosen to be T independent, so that 
the only relevant energy scales are the much lower T c 
and scattering rate. 

For T < T c , inelastic scattering, which has been ig- 
nored in the above calculations, may affect any direct 
comparison between theory and experiment. This leads 
to non-monotonic us dependences and is expected to be 
most important in the optimal and over-doped regime, 
where inelastic effects are strongly T dependent below 
T c , reflecting the T dependence of A(T). It follows, sim- 
ilarly, that in the underdoped regime where A(T) is con- 
stant, inelastic effects should not contribute appreciably 
to AN. 

Finally, it can be noted that the only, nevertheless 
ambiguousc-3, evidence .-for a pseudogap appears in the 
experimental literatures via Drude fits to <j(us), where 
the resulting 1/t(u>) is suppressed at frequency us w 600 
cm . Because of its x independence we, speculate that 
us a may reflect a phonon or other bosonEJ with coupling 
modulated by the opening of an excitation gap. Although 
these secondary pseudogap effects are beyond the scope 
of this paper, we find here that Drude fits are quite gen- 
erally problematic for cases (ii) and (iii) discussed here. 

In summary, in this paper (i) we have investigated 
competing proposals for the pseudogap origin and have 
presented theoretical (and rather generic) mean field ex- 
pressions for o~(us). (ii) Both major mean field pseudo- 



gap schools lead to the prediction that pseudogap ef- 
fects must enter as a rearrangement of excess spectral 
weight (relative to strict BCS theory predictions) in the 
ac conductivity. For the intrinsic school, there is a rapid 
growth in o\ with increasing low T, at fixed low frequen- 
cies (which may have been observed experimentallyci), 
as well as long high us tails, (possibly also observedau) 
which tails become more pronounced with increasing T. 
While for the extrinsic school, there is a transfer of spec- 
tral weight from low to high us. However, here the re- 
sulting secondary maximum is present at all T <T* and 
does not participate in the condensation. Direct com- 
parison with experimental data is made complex by the 
small integrated spectral weight in the condensate and 
our neglect of inelastic scatterings. Therefore, we have 
focused on the condensing portion of the spectral weight 
by considering the change in cri(u>) from to T c . (iv) 
We find that when A(T) w const, strict BCS theory can- 
not account for the relatively high (~ 4A(0)) frequency 
scale associated with the condensate. The extrinsic ap- 
proach suffers from a similar problem. In an artificial way 
rescaled BCS theory remedies this weakness. However, 
the observed high frequency contributions arise naturally 
in an intrinsic model; their presence reflects the under- 
lying scattering of non-condensed bosons from the two 
particle fermionic states. 
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